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I. INTRODUCTION 



In recent years, alternative field theories of gravity (scalar-tensor theories) have been 
analyzed in many contexts (see |1| for a review), specially as a way to study possible 
deviations from general relativity and trough comparison with experimental results put 
bounds on the magnitude of such deviations. 

The basic assumption of these theories is the existence of a single or multiple funda- 
mental scalar fields which are perhaps the effective relics of a fundamental unified theory 
like superstrings, supergravity or Kaluza-Klein theory, and that might couple to gravity 
either minimally or non-minimally. 

One of the most interesting effects occurring in a certain class of scalar-tensor theories 
of gravity is the phenomenon of spontaneous scalarization in neutron stars that was re- 
cently discovered by Damour and Esposito-Farese 0,^, which consists in the appearance 
of a non-trivial configuration of a scalar field in the absence of sources and with vanishing 
asymptotic value. This phenomenon requires the presence of a scalar field nonminimally 
coupled to gravity. This nonminimal coupling can be thought as a field- dependent gravita- 
tional constant in the spirit of the Brans-Dicke theory 0-0]. The occurrence of this highly 
nonlinear effect can be used to put very stringent limits on the deviations from general 
relativity, and in fact using the binary pulsar, the limits for certain range of parameters 
seem to be more severe than those that can be inferred from solar system experiments 
i- 

One puzzling aspect of the phenomenon is that occurs exactly in the regime where the 
Newtonian analysis indicates it should not, while on the other hand, no such phenomenon 
has been said to occur in the regime where the Newtonian analysis indicates it should. We 
give a very short description of these facts based on an energetic analysis and then carry 
out a numerical study to examine its validity and explain why does the full fledge theory 
behave in the opposite direction as compared to what is indicating by the Newtonian 
intuition. To do this we have to perform the analysis in the physical frame where a 
counterpart to the intuitive Newtonian ideas can be found in contrast with the original 
analysis @,@|, which was carried out in a conformal frame that, although yield simpler 
equations, corresponds to a description that is less susceptible of comparison with the 
Newtonian analysis. 

Finally, and in view of the fact that the phenomenon when applied to the binary pulsar 
is being used to put severe limits on the theory's parameters, we investigate the degree 
to which the phenomenon depends on the equation of state assumed for the description 
of the nuclear matter of neutron stars. 

The paper is organized as follows. In section II we briefly review the phenomenon 
of spontaneous scalarization. In section III we provide the naive Newtonian energetic 
analysis of spontaneous scalarization. In section IV we introduce the scalar-tensor theory 
we will analyze in detail, derive the equations of motion for a spherically symmetric, static 
space-time, and make some comments about the realistic equations of state used in this 
work. In section V we discuss the boundary conditions for the numerical integration, 
and define the global quantities to be computed and which are of astrophysical interest. 
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In section VI we provide the numerical analysis and the explanation of the origin of 
spontaneous scalarization on relativistic energetic grounds. Finally, in section VII we 
discuss and summarize our results. 



II. SPONTANEOUS SCALARIZATION 



As we mentioned, the phenomenon of spontaneous scalarization was discovered by 
Damour and Esposito-Farese in scalar-tensor theories of gravity 0,0, and further studied 
in the context of stability analysis [^jTHj- However, we have to comment that Zagaluer 



P] had previously argued for large deviations from general relativity in neutron stars in 
a particular class of scalar-tensor theories. Such results have been criticized in that seem 
spurious due to an artificial definition of some parameters [Q]. 

In these theories, gravity is described both by the standard type (0,2) tensor field, and 
by one or more scalar fields. We concentrate for simplicity in the case of one scalar field, 
and the Lagrangian is then taken to be: 

Where Go is the Newton's gravitational constant, ip denotes the matter fields, and g^i, 
stands for the physical metric (Jordan metric). Using the conformal transformation to a 
non-physical metric gf*^ = f{4>)gfiu (Einstein metric) and redefining a new scalar field. 



9/(0),2^l5(0)^'/' 



4/2(0)^ 50 ^ 2/(0) 
this Lagrangian can be brought in the form 



(2) 



^ = (l^) v^[^* - 2(V.^)'] + C^A^, f-\^)9;.) . (3) 

Here the tensorial operations are performed by using the non-physical metric g*^^. The 
reason given in for writing the field equations in terms of the non-physical metric 
rather than the physical one is the mathematical simplicity related to the fact that the 
non-physical gravitational variables, g*^^, and y?, are minimally coupled. Many authors 
have also employed that frame to analyze neutron-star equilibrium configurations in other 
related contexts PHT2|. 

The field equations derived of the above Lagrangian are 

G;, = 2V;^V> - glA'^.^f + 87rGo/-^(v^)T^. , (4) 
□V = 27rGor'(51n//9<^)T, (5) 



(6) 



where T'^'^ is the physical stress-energy tensor. For a particular class of scalar-tensor theo- 
ries and using this frame (Einstein frame), the effects on neutron stars were investigated. 
Specifically, it was analyzed the equilibrium configurations of a spherically symmetric, 
static, isolated neutron stars using polytropic equations of state to model the nuclear 
matter [Pip|-[T^, and also in axisymmetric slow rotating polytropes P]. 

In particular, it was found that for the case of scalar-tensor theories with l3{(po = 0) < 
—4 [where Pi^p) = {—l/2){d'^\nf/d(f^)], there existed a certain critical baryonic mass 
above which the star developed a non-zero scalar field (p, even if (/jq = (here (po is the 
value at spatial infinity of the scalar field, i.e., the cosmological value) [0,^0. This occurs 
because such configurations are energetically more favorable than the one having ip = 
in all the space (both configurations having equal and constant baryonic masses). This 
striking effect was baptized as spontaneous scalarization in analogy with the spontaneous 
magnetization arising in ferromagnets below the Curie temperature. These equilibrium 
configurations lead to strong deviations from general relativity in the dynamics of binary- 
pulsar systems 0. The fact that the analysis was done in the non-physical frame does 
not allow us to understand the interplay of the true contributions to the energy due to 
the physical fields that are responsible of this phenomenon. We refer the reader to Ref. 
||T3| for a discussion of the problems associated with this point. 



III. NEWTONIAN ENERGETIC ANALYSIS 



One would like to understand the occurrence of the phenomenon of spontaneous scalar- 
ization in neutron stars on simple energetic grounds, much in the same way that one 
understands the occurrence of spontaneous magnetization on ferromagnetic materials by 
comparing the free energy of the states with no magnetization to the free energy of the 
states with magnetization arriving to the conclusion that for a sufficiently low temperature 
the free energy is minimized by the latter rather than the former. 

Let's try then to understand the phenomenon on these energetic terms using a New- 
tonian analysis of the contributions to the total energy of a given matter-scalar field 
configuration. First of all we recall that one expects the true equilibrium configurations 
in this case to correspond to the minima of the energy within the space of configurations 
with fixed values of the other conserved quantities [0, in this case the only relevant 
conserved quantity is the total baryon number of the neutron star (for simplicity we are 
taking the temperature to vanish). From this point of view one would think naively that 
since the net contribution of the scalar field to the energy is expected to be positive, the 
minimizing configuration should have = 0. A little thought then points to a caveat given 
the fact that in the class of scalar-tensor theory with /(0) = (1 + 167r^0^), the effective 
gravitational constant depends on 0, i.e., Ges = Go/{^ + ISvr^^^). The point is that we 
expect the total energy of the configuration to be given by the sum of the mass energy of 
the neutrons (neglecting any thermal energy and those contributions arising from short 
range interactions), the energy of the scalar field, and the gravitational binding energy. 
Thus roughly we expect 
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M = l pU-3^)d'x + I p^ix)d'x - Geff(0) / + ... , (7) 

where we took into account the dependence of the effective gravitational constant on 
0. Now we recall that we are expected to minimize subject to the condition 

A^bar = J Pha.vix)d^X = COUStaut . (8) 

Thus the first term in Eq. (|^ is fixed, and the second is expected to be positive, 
but it is conceivably that if the third integral is sufficiently large (i.e., if the object is 
compact enough) then the positive energy contribution of the scalar field can be more 
than compensated by the change in the Ges that would result in the gravitational binding 
energy becoming substantially more negative, and thus making the configurations with a 
non-vanishing scalar field energetically preferred. This, of course, would require that Gefi 
be an increasing function of 0, which in our specific case corresponds to ^ < 0. Thus we 
seem to have arrived to the conclusion that the phenomenon of spontaneous scalarization 
is possible only if Ges be an increasing function of 0. This is precisely the opposite of 
what is found in the relativistic analysis carried out in Refs. and in the present 

paper. 

We will examine next the phenomenon using the exact form of the theory in order to 
understand why did the naive Newtonian analysis fail so miserably. 



IV. FORMULATION OF THE MODEL 

We will consider a model of a scalar field non-minimally coupled (NMC) to gravity. 
One of the simplest models of this kind is obtained by considering the Lagrangian 



levrGo 



+ c 



mat 



(9) 



Here ^ stands for the NMC constant, and V{(j)) is a scalar potential. In this model the 
schematic matter Lagrangian £mat is intended to represent the ordinary matter contribu- 
tion to the total Lagrangian. 

As we mentioned, equation (|) shows that the introduction of the coupling term is 
equivalent to consider an effective gravitational constant which explicitly depends on the 
scalar field. We note that for ^ positive, (the convention that we keep in all the paper), 
the effective gravitational constant might only decrease with respect to the Newtonian 
value: 



G. 



Gq 



eS 



(10) 



The gravitational field equations following from the Lagrangian (^ can be written as 

1 



where 

rpflU 

Go 



= ^ (4eTr + + TZt) , (12) 



Tf = V'^(0V'^0) - (7^'^Va(0V^0) , (13) 
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(14) 



The energy-momentum tensor of matter T^^^t will be represented by a perfect fluid de- 
scribing the cold catalyzed matter of a neutron star ||15|| : 



TZt = {p + p)U^U''+pg^'' , (15) 

which possesses the symmetries of the spacetime. The scalar field will also be assumed to 
posses these symmetries. 

Finally, the equation of motion for the scalar field becomes 

□0 + 2?^fi=^. (16) 

We will focus on a metric describing spherical and static space-times. Four our pur- 
poses it will be convenient to adopt the so-called radial gauge coordinates which are also 
maximal slicing coordinates (hereafter we refer these coordinates as to RGMS coordi- 
nates) : 

ds^ = -N^{r)df + A\r)dr^ + r^dO^ + sin^ Od^^ . (17) 

Our goal is to analyze solutions of the gravitational, matter, and scalar field equations 
describing neutron-star models and the resulting space-time. Owing to the complexity of 
the resulting equations, we will perform a numerical analysis. To do so it will be helpful 
to adopt the following variables 

u{r) = ln[A^(r)] , (18) 
z>(r) = z/(r) - z/(0) , (19) 

^(,)=(l-H^)"'^ . (20) 
The relevant Einstein equations take then the following form 

= , (21) 

| = A^{^ + 4.rG„ri,} . (22) 

where 

E = N'^Tl^ , (23) 



is the effective total energy density. 

On the other hand, the Klein-Gordon equation can be written directly in terms of 
sources as follows: 



□0 = -levr^^Go {E~S) + 



dV{(l)) 



(24) 



where 



S = T, 



eff i 1 



(25) 



is the trace of the spatial part of T^f^ , which plays the role of an effective pressure. 
In the RGMS coordinates this equation reads 



- + — - ( 1 

r dr 



+ 1 



1 



} 



2Gom 



167r^0Go (E - S) 



90 
dr 



The matter equations are obtained from 







(26) 



(27) 



and straightforward manipulations with help of the Einstein equations show that the 
matter and the scalar-field energy-momentum tensors are conserved separately, leading 
thus to 



Y7 r-piiv 
^ puJ- matt 







(28) 



This yields the equation of hydrostatic equilibrium which is similar to the one obtained 
in general relativity with no scalar fields. 

For cold catalyzed matter (see below), it's customary to have the equations of state 
(EOS) parameterized by the baryon density n(r), thus it is convenient to write the equa- 
tion of hydrostatic equilibrium in terms of this: 



dn 
dr 



(p -|- p)n dv 
7P dr 



(29) 



^llj^ is the adiabatic index. 



where 7 : = 

It is clear from Eqs. (p!2| -|l4D, that the intermediary variables E and S [see Eqs.(p3|), 
(p5|)] involve second order derivatives of the scalar field. However, we can eliminate such a 
terms from the gravitational field equations with the help of the Klein-Gordon equation, 
and obtain "sources" containing at most first order derivatives of the scalar field. We also 
introduce the following dimensionless quantities (where we restore the factors of c) 
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m 



n 



m 

Mo 
n 

.Ifm 
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p := p ■ 



'Gq 



"'^ ■ 1 + 1671^02 ' 

then the final form of the equations to be solved numerically is : 



— + 47rrGefr 



1 + 1677^^0 ((9,.0)Geff 

'^ + dfD-(l-—Y' [AnfE- 
r \ f 



{df(j)) -V{(t))+p 



2^2 

m 



(30) 
(31) 
(32) 

(33) 
(34) 
(35) 
(36) 
(37) 

(38) 
(39) 



^ 2m Y 

f / 



dV{<p) 



IQ-Kicf) (e - S) 



dffi = 
where 
E-S 

E 



{p + p)n 



IP 



(40) 
(41) 



off 



1 + 19277^202^-^^ 

Gcs 

1 + 19277^202^-^^ _ 

dV{(f)) 



A^ 



{df4>)\l + 120 + 4V^(0) +P-Sp + 12^0 



(42) 



A2 



(l + 1927reV'G'efr) + ^{dr<P? (1 + 8^ + Q47ie<P'G,s) 



+4e</'^|^ + ^(0) (l - 6471^^'^^) + p (l + 1287reV'Geff) + 1927ie^^G,sP 

Here and 5* are dimensionless as in Eq. (|33|) , and we have dropped-out the tilde over 
0. 

We note that (^) with (^) generalizes the Volkoff-Oppenheimer equation of hydro- 
static equilibrium for the model (|^). We also note that the sources of the differential 
equations contain only first order derivatives of the field variables and are thus suitable 
for numerical integration with a Runge-Kutta algorithm. 



(43) 



By subtracting Eq. (^) from (^31) we can obtain an effective total pressure S/3 from 
wliich we recognize individual contributions, like the pressure terms associated with the 
scalar field 

1 ,(9,0)2-^(0) , (44) 



and the NMC contributions 

Gcs 



4 



3(1 + 1927r^202(5^g) 
1 



6^2 



2^2 



(5,0)^ 



1 - 16^ + 64<V^G'cfr 



1 + 1927r^V^G'cfr 
dV{(j)) 



y(0) (3 + 647reV'G'eff 



Pmat 



G. 



cS 



3(1 + 19271^202^'^^^ 



128n^^ (p^G^sp + 3p(l + 647r^V 



For simplicity we will consider hereafter only the case V{(f)) = although the basic 
results are expected to remain qualitatively unchanged. 

Regarding the equation of state, we make the standard assumption that the nuclear 
matter in neutron stars is considered to be at zero temperature (cold catalyzed matter) 
| |15| |. This matter is represented by an equation of state usually parameterized by the 
baryon density number: 



p = p{n) , 
P = p{n) ■ 



(47) 
(48) 



In general the EOS are not given analytically. In particular, for realistic EOS, i.e., those 
built from effective field-theories or many-body nuclear calculations, the above set of 
equations is given in tabulated form. Thus, a preliminary interpolation is required. In 
this paper we will use the collection of EOS of |T6|JT7[] constructed with a logarithmic 
interpolation. In particular we show results for three EOS representing the different de- 
grees of stiffness typically found in other models. These are the model of Pandharipande 
in which only neutrons are taken into account [|19| (hereafter referred as to PandN; rep- 
resentative of a "soft" EOS), the model II of Diaz-Alonso 0] (hereafter referred as to 
DiazII; representative of a "medium" EOS), and the model "0.17" of Haensel et al. 
(herafter referred as to HKP; representative of a "stiff" EOS). We point out that the 
DiazII and HKP models, being relativistic models, are thus causal in all the domain of 
baryon densities. The PandN model, despite being non-relativistic, is causal for the range 
of densities usually found in stable configurations of neutron stars. We refer the reader 
to the previous references, and the bibliography therein for a more detailed description of 
such EOS. 



, (45) 
(46) 



V. BOUNDARY CONDITIONS AND NUMERICAL METHODOLOGY 



Interior solution. The regularity condition at r = (center of the star) on the 
metric requires the boundary condition 



Q 



m(0) = 



(49) 



The boundary condition on z/(r) is by definition [cf.Eq. (19)] 

i>(0) = . (50) 

The integration of Eq. (^) is performed by specifying the baryon density at the center 
of the star 

n(0)=ne , (51) 

with the regularity condition drn{0) = 0. 

In a similar way the boundary and regularity conditions on the scalar field are 

m = <pc , (52) 

dr(f){0) = . (53) 

On the other hand, the value (pc cannot be arbitrary, but must be so that satisfy the 
appropriate boundary conditions at spatial infinite. This is enforced by the use of a 
standard shooting method |^l|]. In the present paper, we only consider an asymptotically 
vanishing scalar field. This corresponds to the situation in which spontaneous scalarization 
might arise. The situation for a non-zero asymptotic value of the scalar field can give rise 
to "induced" scalarization in the neutron star, and can be also obtained in exactly the same 
fashion by a shooting method. Figure |l] shows the solution of the Klein-Gordon equation 
satisfying vanishing asymptotic boundary conditions. Figure ^ shows the "trajectories" 
(0c, Pc) for which those conditions are verified within different neutron star configurations. 

Exterior solution. We can compute the exterior solution by integrating the equa- 
tions from the star surface to spatial infinite. Outside the support of the fiuid variables p, p, 
we "compactify" the space by a transformation on the r coordinate of the kind u = 1/r, 
and integrate the resulting equations (see the Appendix) in the domain u G [1/i?, 0]; at 
the star radius R (which is obtained numerically) and which is defined as the coordinate 
r at which the pressure vanishes p{R) = 0, the values of the field variables are given 
numerically, and they represent the boundary conditions for the integration of the equa- 
tions outside the star. The integration constants are fixed by matching continuously the 
interior with the (numerical calculated) exterior solution. 

The physical lapse N{r) = e'^ at r = is calculated at the end of the numerical 
integration from (p!9|) 

^^(0) = -z>oo , (54) 

where the value i>oo is obtained from the numerical integration, and it ensures that at 
spatial infinite our coordinates correspond to the standard Minkowski coordinates. Figure 
^ depicts the metric potentials with asymptotic conditions matching with Minkowski 
coordinates at spatial infinite. 

Global quantities. The staticity of the configurations we are studying ensures that 
the ADM (gravitational) mass and the Komar mass coincide, and can be easily evaluated 
from the integral 

in 



M, 



ADM 



POO 

liiiir-^oo 'm{r) = An r^E{r)dr . (55) 

JO 



We emphasize that the energy density E might not have compact support as it might 
include contributions of the scalar field. Therefore in order to compute the actual ADM 
mass the integration have to be performed from the center of the star to spatial infinite. 

Moreover, the conservation of the baryon number leads to the conserved total baryon 
number of the star given by 

J\f = Atx nir)Air)r^dr . (56) 
Jo 

The total baryon mass is defined by 

Mbar = mfcAr , (57) 

where irih is the mean mass of single baryons (1.66 xlO^^'^ kg). 

As it is usual, we can define the total binding energy of the star as 

^bind = ^ADM — ^bar , (58) 

and the fractional binding energy as 

^^frac=l-^ . (59) 

From the asymptotic behavior of the scalar field with vanishing boundary conditions, 

<P{r)^^ + 0{l/r') , (60) 
c r 



one can define the so called "scalar charge" uj as follows: 

^2 



uj := — lim 





Go dr 



(61) 



Following 0,0, we introduce the coupling strength 



JWadm 

In the next section we provide the numerical analysis of the above quantities, and 
their dependence on the parameters pc = pi^c), (f>c, and the equation of state. 
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VI. NUMERICAL RESULTS 



A. Relativistic energetic analysis of spontaneous scalarization 



In section III we showed that a naive Newtonian analysis leads to the conclusion that 
the theory described by with ^ > should not give rise to a spontaneous scalarization 
since the presence of the scalar field would tend to increase the total energy by decreasing 
of the negative binding energy and thus to a configuration with a larger total energy 
than the corresponding case in absence of scalar field, for a fixed baryonic mass. In this 
section we will show that a numerical analysis of the full relativistic energy contributions 
confirms that the phenomenon of spontaneous scalarization indeed occurs and that the 
naive analysis fails. 

Let us separate the total energy density E [see Eq. (^31)] into contributions that can 
be identified as representing the different effects that arise in the theory. This separation 
of course can be consider to have no more than heuristic value, for even the energy 
density itself has no invariant meaning in a diffeomorphism invariant theory. We do this 
separation in the following way 



Pint l?^ J 



El-- 



p — mi,n{r) , 
mhn{r)A{r) 



2^2 



G, 



eflf 



1 + 1927r^202^^^ 



(l + 1927r^V^G'cff 



— Ej-^ 



G. 



eflf 



1 + 192^202^^^ 

Then it is clear that 



p (l + 1287reV'Gefr) + 192tt ^fl^^G^sP 



(63) 
(64) 
(65) 

(66) 



(67) 
(68) 



E 



Pbar 



PhZ'^ + Pint + E^ 



El 



(69) 



This decomposition has the following interpretation: p^^j. is the physical baryon density, 
that is, the energy density that will give rise to the total baryon mass; pmt is the inter- 
action energy between baryons (usually it is given in a complicated way related to the 
physics model of the nuclear matter), — Pb™*^ is the negative standard binding energy- 
density of baryons (is the part that is associated only with the difference of the baryon 
density multiplied by the geometrical factor corresponding to the proper volume element 
on spacelike hypersurfaces and the baryon density itself. The integral of this difference 
is what yields the total binding baryon energy in the case of pure general relativity |p22||), 
Eff, is the energy- density of the scalar field as it is usually defined without NMC, E^ is the 
contribution to the energy-density of the scalar field due to the non-minimimal coupling 
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(note that this vanishes when ^ = 0), and finally, is the contribution to the total 
matter-energy density due to the presence of the non-minimal coupling (we also note that 
it vanishes for ^ = 0). The integration of each of the terms from r = to r = +oo will 
give the total energy contributions: 

MadM = Mbar + Mint - Mbind + + M| + • (70) 

We remark that the above interpretation is sustainable as far as we concern ourselves with 
the functional form of the terms themselves and regard them as applicable to different 
test configurations. On the other hand the particular configuration that minimizes the 
total Madm arises from an interplay of all such effects and therefore the actual value of 
each of the terms in Eq. (|70D is affected by the presence of the other terms. Nevertheless 
the heuristic value of the arguments is validated since such affectations can be consider 
as higher order effects. 

In the absence of scalar field we have 

M^DM = <a^ + Mnt''-<nd ■ (71) 

So if we want to compare two configurations at fixed baryon mass (Mbar = M^^), one 
with 7^ and the other with null scalar field, we have 

AM := Madm - M^^^ 

= Mint - M:1^ - (Mbind - M^^^) +M^ + Ml + Ml . (72) 

Now according to the Newtonian intuition we expect that |M^^| > |Mbind|- So the term 
— ^Mbind — M^^^^ would be positive. Moreover the energy M^ is also positive and we 
expect that M^t — is small. In order to have spontaneous scalarization we need 

AM < and so the term M| + M| should be negative enough to more than compensate 
for the positive energy contributions and also for the decrease in the ordinary binding 
energy. Figure ^ shows the local behavior of such energy contributions (referred with 
lower-case letters in the caption) for a configuration with spontaneous scalarization, and 
confirms the previous conjectures. Figure |^ depicts the corresponding energy-densities 
as defined by the Eqs. (|63|)-(|68D. Table I shows examples of configurations exhibiting 
spontaneous scalarization as compared with those in absence of scalar field at the same 
baryon mass. We can appreciate from table I that the different energy contributions add 
to reproduce the ADM mass and that the value for the configuration with 7^ is lower 
than its corresponding value with = 0. Note also that the binding baryon energies 
in absence of scalar field are smaller than their corresponding in the configurations with 
(f) ^ 0. This confirms the original Newtonian intuition. However, the relativistic energetic 
analysis also shows that the Newtonian naive argument failed to explain the spontaneous 
scalarization with ,^ > since it did not take into account the large negative energy 
contribution M| which is the responsible for the decrease of the total energy (M| turns 
to be always positive), and as it is shown in the examples of table I, is the dominant 
part of the above difference AM. This contribution as well as M| have no Newtonian 
counterpart. 

1 



B. Effects of the equations of state 



A simple and enlightening way to present the effects of the equation of state on the 
phenomena behavior of the other global quantities as a function of the total baryon mass 
of the neutron star. Figures ^ and ^ show the behavior of the coupling strength a [see 
Eq. ( |B^ as a function of the total baryon mass for three EOS. We note that there are 
critical values beyond which the phenomenon of spontaneous scalarization ensues. 

The critical values correspond to the point at which configurations with a non-vanishing 
scalar field become energetically more favorable than the configurations with 0(r) = 0. 
This behavior can be better appreciated in Fig. |[ It is interesting to note that M^"* 
is an increasing function of the stiffness of the EOS (for example when = 2, M^"j.* for 
the HKP EOS increases in ~ 37% relative to that for the PandN EOS.), but that the 
maximum values of a seem to be almost independent of the EOS. 

Figure |^ shows that the larger the value of ^ the less important the effect of the EOS 
on the critical value M^^. This can be understood on basic grounds by looking at the 
expression of E^^ [Eq. (^)], which as we discussed is the main responsible for spontaneous 
scalarization. As ^ increases i?| —p. This means that its relative contribution is be- 
coming dominant regardless of the central- value of p. Therefore spontaneous scalarization 
can take place for smaller values of pc- Furthermore, the lower the value p^ the more alike 
are the equations of state since the three EOS are built with the same standard neutron- 



star-crust EOS which corresponds to that of So the configurations of the star with 



lower central values of p are very much alike. This is the reason for the critical baryon 
masses look very similar for different EOS with large values of ^. Table || shows the values 
of quantities corresponding to the onset of spontaneous scalarization (i.e., corresponding 
to the critical Mbar)- 

As it was pointed out in |p, the non-perturbative effects of the scalar field in the 
neutron star are the analogue of spontaneous magnetization in ferromagnets at temper- 
atures lower than the Curie point. In the Landau ansatz where the free energy F of the 
ferromagnet is a function of the magnetization p in the form F = a(T — Tc)p? j2 + bp'^/A, 
it appears two non-trivial minima for temperatures T < Tc- The analogue to the curve 
giving the locus of the minima of F, is depicted for the case of spontaneous scalarization 
in figure |^, where the scalar charge u plays the role of p, the total mass Madm the role of 
F, and finally the baryon mass is the analogue of the temperature. Obviously in this case, 
the complexity of the system makes difficult an attempt to provide an analytic expression 
of the form Madm = Madm{^, Mbar) for which we could recover the behavior of figure ||. 

Another effect related to the phenomenon of spontaneous scalarization is the increase 
of the maximum masses (both Madm and Mbar) as compared with the values allowed by 
general relativity. Figure |10] shows the dependence of Madm on p^ for the selected EOS. 
The behavior of these profiles are qualitatively the same as in general relativity (see Ref. 
|T6|). Table |ITl| gives the relevant properties of the neutron star models corresponding to 
the maximum mass configurations. Note the increase of Madm with ^. 
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VII. DISCUSSION 



GM = G'eff(0) 



One of the main goals of this work was to understand the discrepancy between the 
Newtonian expectations and the relativistic energetic analysis. A fortiori it turns out 
that the difference can be understood by noting the different roles that the gravitational 
constant plays in the energy considerations in the Newtonian and the relativistic theories. 
As discussed in section 111, the Newtonian point of view leads to the expectation that the 
total mass should be expressed in the form: 

M^J p^U^)d'x + J p^{x)d'x - CM) J ^^^^^^^^d'xd'y + ... , (73) 

so that the gravitational constant appears as the coefficient of a term that represents the 
binding energy which is negative, so a decrease in Ges results in an increase on Ad. In the 
relativistic theory the ADM mass itself is of gravitational origin because it is a measure 
of the departure of the metric (at large distances) from the Minkowski metric, i.e., it 
represents the 1/r coefficient of such deviation. In this sense we can say that the only 
thing which has meaning is the combination G x Ai and not Ai itself. Thus the quantity 
that must be minimized is (heuristically) : 

\x-y\ 

(74) 

Which we can write as 

GM = Geff (0)A^ba. + Geff (0)A1<^ - GI^ I ^^^^^j=^^^dVy + ... , (75) 

J \x-y\ 

The point is that the appearance of a nonzero value of (f) which results in the second 
term becoming positive (from a zero value when = 0) can be more than compensated 
by the decrease in the first term due to the decrease in Ges that is associated with an 
increase in 0. The change in the third term, which in effect is in the direction inferred from 
the Newtonian analysis, playing now a secondary role. This explains why the Newtonian 
analysis failed and why the relativistic analysis shows that the phenomenon of spontaneous 
scalarization is associated with an effective gravitational constant that decreases with 
0. This explanation is supported by the numerical analysis since the term M| which 
represents precisely the effect described here turns out to be the one that is responsible 
for the lowering of the ADM mass of the configuration with a scalar field in comparison 
to the configuration without one. 

The second goal of this paper was to study the dependence of the parameters associated 
with the spontaneous scalarization with the equation of state of the nuclear matter. This 
can be summarized as follows. First the critical baryon masses at which spontaneous 
scalarization develops increases as the equation of state becomes more rigid. As ^ increases 
those critical values become much less dependent of the EOS. Second, the maximum values 
of the coupling strength a seem to be independent of the EOS. Finally, the effect of the 
scalar field makes the neutron star to support larger masses than in general relativity, 
and this is more evident as ^ increases. 
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APPENDIX: DIMENSIONLESS FORM OF FIELD EQUATIONS AT THE 

EXTERIOR 



By performing a "conformal" transformation of the kind m = 1/f in Eqs. (p^)-P3|) 
(with p = p = V{(f)) = 0) we obtain the equations to be solved in the domain u E [1/R,0] 
which corresponds to the range from the radius of the star to spatial infinite. 



duTn 



u 1 
1 - 



2mu) 



-1 



771 + 47rGcfr(l — 2m-u) 



1 - 16TT^U(f){du(l))GcS 

— (1 — (m — AiruEy^ du 



u 



1 + 1927r^202^^^ 

1 + 19277^202(5^^ 

+ ^(9„0)2(l + 8e + 64<VG'efr 



(9„0)2(1 + 12O , 

-A^(t>{du(i>){duV) (l + 19277^^^5 



(Al) 
(A2) 

(A3) 



(A4) 
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TABLES 

TABLE L Examples of configurations exliibiting spontaneous scalarization compared with 
those at the same baryon mass in absence of scalar field. is the non-minimal coupling con- 
stant; Mbar is the total baryon mass; Mbin is the baryon-binding energy; Mint is the total 
baryon-interaction energy; is the total mass contribution due to the single scalar field; M| 
is the mass contribution of the correction to the total matter-energy by the presence of the 
non-minimal coupling; m| is the mass contribution of the correction to by the presence of 
the non-minimal coupling; Madm is the ADM-mass. We use Mq = 1.989 x 10^*^ kg. 



EOS 




4>{r) 


Mbar 


Mbind 


Mint 


M^ 


M« 


m| 


Madm 








[Mo] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


[Mq] 


DiazII 


2 




1.563 


-0.179 


0.082 


0.011 


-0.226 


0.179 


1.429 




2 


= 


1.563 


-0.207 


0.077 











1.433 




6 




2.170 


-0.320 


0.170 


0.011 


-0.559 


0.413 


1.885 




6 


= 


2.170 


-0.469 


0.204 











1.904 



TABLE IL Critical configurations towards spontaneous scalarization. ric, pc are the central 
baryon and proper energy densities respectively, Pc is the central pressure, R is the radius of the 
star, Nc is the value of the lapse at the center of the star. The remaining symbols are defined 
in the caption of Table |. We use Go = 6.672 x 10"" m^kg-^s'^ c = 2.9979 x 10^ ms-^ 
Mq = 1.989 X 10^° kg, and pnuc = 1-66 x lO^'^kgm"^ 



EOS 




ric 
[.lfm-3] 


Pc 
[PnucC^] 


Pc 
[PnucC^] 


Madm 

[Mq] 


Mbar 

[Mq] 


R 

[km] 




SGqMadm 


Nc 


-'W'adm 


c-'R 


PandN 


2 
6 


4.22 
1.76 


4.46 
1.80 


0.254 
0.026 


0.596 
0.1779 


0.617 
0.1779 


11.29 
18.65 


-0.035 
-5.87 xlO"^ 


0.156 
0.028 


0.822 
0.944 


DiazII 


2 
6 


2.52 
1.17 


2.65 
1.20 


0.150 
0.018 


0.757 
0.2077 


0.787 
0.2079 


13.90 
20.08 


-0.038 
-0.0012 


0.160 
0.030 


0.824 
0.945 


HKP 


2 
6 


1.97 
1.243 


2.05 
1.27 


0.122 
0.020 


0.811 
0.1970 


0.845 
0.1973 


14.24 
18.98 


-0.042 
-0.0012 


0.169 
0.030 


0.826 
0.945 
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TABLE in. Maximum mass models. (j)c is the central scalar field, M^. is the "gravitational 
mass" at the surface of the star, u is the scalar charge, a is the coupling strength between the 
scalar field and gravitation. The remaining identifiers are defined in the caption of Tables | and 

The configurations in absence of scalar field corresponds to the maximum mass models in 
general relativity. 



.Ifm 



Pc Pc (t>c Ms Madm 

[PnucC2] [pnucC2] [Mq] [Mq] 



Mbar R iO 1 

[Mq] [km] [Mq] 



2GoMadm 



17.36 24.69 11.69 1.662 1.662 
2 15.51 20.98 8.73 0.0171 1.620 1.694 
6 15.51 20.98 8.73 0.0169 1.570 1.776 



1.932 8.52 -0.162 0.575 
1.981 8.91 0.0536 -0.169 0.558 
2.118 9.13 0.0830 -0.192 0.574 



11.14 15.20 5.30 1.933 1.933 
2 10.75 14.49 4.90 0.0226 1.868 1.990 
6 10.94 14.85 5.10 0.019 1.801 2.085 



2.210 10.92 -0.143 0.522 
2.297 11.21 0.0818 -0.153 0.524 
2.451 11.42 0.115 -0.175 0.538 



6.32 8.76 4.29 2.835 2.835 
2 5.64 7.42 3.27 0.0146 2.768 2.879 
6 5.56 7.28 3.16 0.015 2.667 3.018 



3.421 13.67 -0.206 0.612 
3.492 14.15 0.0805 -0.212 0.600 
3.732 14.56 0.132 -0.236 0.612 



a 



Nr. 



0.310 
-0.0316 0.344 
-0.0467 0.337 

0.375 
-0.0411 0.380 
-0.055 0.368 

0.301 
-0.0279 0.332 
-0.043 0.329 
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FIGURES 



T — I — ' — I — ' — I — I — ' — ' — ' — I — I — ' — ' — I — r 




Q I , , , , 1 , , , , 1 , , 1 ^ r— I , , 1 , , 1 

12 3 4 

logio[l + r/km] 

FIG. 1. Scalar-field as a function of r for = 6, ric ~ 0.7 fm~^ (pc ~ 8.37 pnucC^) and 
(pc 0.0221 for the EOS DiazII. The asterisk depicts the location of the star surface {R ~ 12.63 
km). The computations includes the domain from r = to r = +oo, however for pictorial 
convenience here the logarithmic scale shows the plot only to some extent which includes various 
orders of magnitude beyond the star surface. 



on 



5 10 15 20 25 30 35 40 

P [P c^l 
r'o Lfnuo J 

FIG. 2. Functional dependence of (pc vs pc obtained from a shooting method with = 2 for 
three different equations of state: HKP (dash-doted hne), DiazII (sohd hne) and PandN (dotted 
Une). The appearence of the non-trivial solutions {(pc / 0) marks the onset of spontaneous 
scalarization. 
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FIG. 3. Metric potentials as a function of r corresponding to the case of Fig. |. The curves 
refer to the lapse function N{r) (solid line), to the metric potential A{r) defined by Eq. ([20| ) 
(dashed line), and to the product AN (dash-dotted line). Note that unlike the Schwarzschild 
exterior solution this product differs from unit outside the star. Here the asterisks indicates the 
location of the star surface. 
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CO 




logioLl + r/km] 

FIG. 4. Energies for the same case of Fig. |l| as a function of r. The profiles represent 
the integration of the energy-densities of Fig. |5| from r = to oo. The upper sohd hues 
(positive values) correspond to m(r) (larger values) and ^^(r) (smaller values) respectively; 
the dashed line stands for m\,s.j-(r); the dotted line corresponds to rrii^tif)', the dash-dotted 
line (positive values) corresponds to m^{r); the dash-dotted line (negative values) stands for 
mbind('')'j the lower solid line (negative values) refers to m^(r). The asymptotic values of these 
profiles correspond to those of table |. Note that the mass m(r) continues to grow outside the 
star (the asterisk depicts the star surface) due to the contribution of the energies associated 
with the scalar field. 
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T 




I I I I I I I I I I I I I I I 

5 10 15 

r/km 

FIG. 5. Energy-densities as a function of r for the same case of Fig. || (see the text for 
description) corresponding to (from the upper to the lower curve respectively) (solid line), 
E (dash-dotted line; larger values), p-mt (dash-dotted line; lower values), (dash-dotted line; 
lowest positive values), p^™'^ (dashed line), and E^ (dotted line). In this scale, the energy-density 
Efj) lies almost on the r-axis. 
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T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.5 1 1.5 2 2.5 3 3.5 4 



FIG. 6. Coupling strength vs the total baryon mass for the EOS labeled HKP (dash-dotted 
line), DiazII (solid line), and PandN (dotted line) for ^ = 2. The dashed segments of the lines 
refer to the unstable configurations. The positive branches {—a > 0) correspond to configura- 
tions with (f> > 0, while the negative ones (—a < 0) refer to those with </> < 0. The critical 
baryon masses for spontaneous scalarization are shown in Table ||. 



OK 



0.5 1 1.5 2 2.5 3 

Mbar [Mo] 

FIG. 8. Total fractional binding-energy as a function of the total baryon mass for the EOS 
DiazII with ^ = 6. The solid line refers to configurations computed with ^ 7^ while the dashed 
line depicts the configurations with ^ = 0. Here second branches (unstable configurations) are 
also shown. Beyond some critical baryon-mass the configurations with 7^ are energetically 
more favorable than the corresponding at the same baryon mass with (p = 0. Note the increase 
in the maximum baryon-mass from 2.21 (cp = 0) to 2.45 {(p ^0). 
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FIG. 9. Scalar charge vs A/adm for the EOS labeled HKP (dash-dotted line), DiazII (solid 
line), and PandN (dotted line) with ^ = 2. The dashed segments of the lines label the unstable 
configurations. 



08 



10 15 20 25 30 35 
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FIG. 10. Total gravitational mass as a function of the central energy-density for the EOS 
HKP (dash-dotted line), DiazII (solid line) and PandN (dotted line) for ^ = 6. The dashed 
segments of the lines correspond to unstable configurations. 
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